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Abstract Credit estimation and bankruptcy prediction methods have utilized 
Altman’s Z-score method for the last several years. It is reported in many 
studies that Z-score is sensitive to changes in accounting figures. Researchers 
have proposed different variations to conventional Z-score that can improve the 
prediction accuracy. In this paper, we develop a new multivariate nonlinear 
model for computing the Z-score. In addition, we develop a new credit risk 
index by fitting a Pearson type 3 distribution to the transformed financial 
ratios. The results of our study have shown that the new Z-score can predict 
the bankruptcy with an accuracy of 98.6% as compared to 93.5% by Altman’s 
Z-score. Also, the discriminate analysis revealed that the new transformed 
financial ratios could predict the bankruptcy probability with an accuracy of 
93.0% as compared to 87.4% using the weights of Altman’s Z-score. 

Keywords credit risk • bankruptcy • prediction • Pearson type 3 distribution • 
Z score • non-linear models • Type II errors • Type I errors 


1 Introduction 

Credit ratings have become an integral part of today’s capital markets as they 
help in the evaluation and assessment of credit risk, benchmark issues and 
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create secondary markets for those aspects. Credit risk exists virtually in all 
income-producing activities and their inappropriate evaluation or inadequate 
mitigation would result in failure of institutions. In general, the credit ratings 
given by agencies such as Standard and Poor, Moody and Fitch are based on 
the probability of default and recovery rate taking into account not only the 
variables in the financial statement of the firms but also the market cues. Pre¬ 
dicting bankruptcy for firms in financial distress from the financial statement 
history is an important problem studied widely by the researchers (Bartual et 
al. 2012; Hernndez & Wilson 2013; Mendes 2014; Zaghdoudih 2013). Among 
them the Altman’s Z-score (Tony et al. 2005; Radu et al. 2009; Altman 1968; 
Altman et al. 1977) is the most popular and widely accepted metric for pre¬ 
dicting the bankruptcy. The popularity of Z-score may be attributed to its 
simplicity in computation and ease in its application (Ali & Kim-Soon 2012; 
Khalid et al. 2008; Allen et al. 2006; Landsman et al. 2009). 

Altman’s Z-score uses mainly accounting figures in the financial statement 
as variables in the computation. The Z-score is highly sensitive to small vari¬ 
ations in these figures due to its dependency on them. This leads to an ex¬ 
aggerated Z-score in case they are manipulated as it does not include the 
past accounting profile of the business into consideration during its compu¬ 
tation. Therefore, the bankruptcy probability predictions using Altman’s Z- 
score would cause significant levels of type-I errors (classifying bankrupt firms 
as non-bankrupt). Beaver et al. (2009) demonstrated the bias of Z-score in 
predicting the bankruptcy. In addition, the models have the weakness of not 
being immune to false accounting practices (Aasen 2013). It is stated by Alt¬ 
man that the retained earnings account is subject to manipulation via corpo¬ 
rate quasi-reorganizations and stock dividend declarations, which may cause 
a bias (Altman 2000). Moreover the weights used by Altman are still preva¬ 
lent even though the financial reporting environment has changed drastically 
from a rule-based approach to a principle-based set of standards aiming for 
harmonization with the International Financial Reporting Standards (IFRSs) 
(Benston et al. 2006; Karim & Tan 2010; Jamal et. al. 2010). 

The Z-score, which is a derivative of financial ratios may not represent 
different risks using same quantitative figures or same business risks for dif¬ 
ferent financial statement figures. To account for this asymmetry, the Z-score 
may be adjusted by including the earnings management in the computation 
procedure (Seong et al. 2012) but it may also be exposed to manipulations 
in the accounting data. These models should not be applied to financial firms 
due to their frequent use of off balance-sheet items (Altman & Edith 2006). 
Also, the results of the model may vary over time, which may be explained 
by the uncertainty of the stock prices as they are subject to the stock market 
opinion. During periods when the stock market is relatively high, the Z-score 
outcomes will be higher than in times when stock prices are low. 

The financial scores are linearly combined to obtain the Z-score using the 
weight functions derived from the multivariate discriminate analysis. In real¬ 
istic scenarios the financial ratios used as independent variables may not be 
linearly related. Also, the score is biased to small variations in the financial 
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scores. Moreover, it is not possible to compare the performance of different 
firms such as non-manufacturing, manufacturing as the weights of the finan¬ 
cial ratios would differ among different firms. Further, it is difficult/not feasible 
to develop specific models tailored to address the scenarios of each type of in¬ 
dustry (retailers, airlines etc.) even though it may look ideal (Altman & Edith 
2006). 

Keeping in view the shortcomings of Altman’s Z-score and the adjusted 
Z-score (Seong et al. 2012), we frame the following objectives for the present 
study: 

1 . developing a score based method using a nonlinear form of financial ratios; 

2 . designing an index using an equi-probability transformation by fitting a 
Pearson type 3 (P3) distribution to the newly developed Z-score say Z M : 

3. formulating a rating scheme based on the index; 

4. comparing the Zm with Altman’s Z-score and the proposed rating scheme 
with those being employed widely by financial institutions for bankruptcy 
predictions. 

The present paper is organized as follows. In Section 2 we propose a new 
nonlinear transformation model for computing the Z-score while in Section 3 
a new index (based on the Zm) using a P3 distribution is developed. The 
methodology followed for predicting the bankruptcy probability of firms is 
presented in Section 4. The datasets and the results are described in Section 5. 
Conclusions and discussion are deferred to Section 6 . 


2 A generalized nonlinear score for modeling the financial ratios 


The standard Z-score is a statistical measurement of a scores’s relationship to 
the mean in a group of scores generally measured by the formula 


Z = 


(X-n) 


(1) 


in which X denotes the set of measurements, /i and a respectively denote the 
mean and standard deviation of the data in the set X. The Z-score is a very 
useful statistic for obtaining the probability of a score occurring within a nor¬ 
mal distribution and comparison of two scores that are from different normal 
distributions. Altman (1968) first proposed a Z-score measuring a company’s 
financial strength using a weighted sum of several factors among the variables 
(financial ratios) that gives an approximate description of the bankruptcy 
probability. 

Altman (1968) utilized a data set composed of sixty-six corporations with 
thirty-three firms in each of the two risk groups and financial ratios given in 
Table 1 to obtain a set of ratios that influence the bankruptcy prediction. The 
mean asset size of these firms was $6.4 million, with a range between $0.7 
million and $25.9 million. Also, the financial ratios R31 (aq), i?26 (# 2 ), -R21 
( 2 : 3 ), .R15 ( 24 ) and 1?19 ( 25 ) were identified as key variables for bankruptcy 
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Table 1 Financial Ratios previously used in bankruptcy prediction studies 


Ratio# 

Name 

Type 

R1 

Cash/Current Liabilities 

liquidity 

R2 

Cash Flow/Current Liabilities 

liquidity 

R3 

Cash Flow/Total Assets 

liquidity 

R4 

Cash Flow/Total Debt 

liquidity 

R5 

Cash/Net Sales 

liquidity 

R6 

Cash/Total Assets 

liquidity 

R7 

Current Assets/Current Liabilities 

liquidity 

R8 

Current Assets/Net Sales 

liquidity 

R9 

Current Assets/Total Assets 

liquidity 

R10 

Current Liabilities/Equity 

liquidity 

Rll 

Equity/Fixed Assets 

solidity 

R12 

Equity/Net Sales 

solidity 

R13 

Inventory/Net Sales 

liquidity 

R14 

Long Term Debt/Equity 

solidity 

R15 

Market Value of Equity/Book Value of Debt 

solidity 

R16 

Total Debt/Equity 

solidity 

R17 

Net Income/Total Assets 

profitability 

R18 

Net Quick Assets/Inventory 

liquidity 

R19 

Net Sales/Total Assets 

profitability 

R20 

Operating Income/Total Assets 

profitability 

R21 

Earnings Before Interest & Tax/Total Interest Payments 

liquidity 

R22 

Quick Assets/Current Liabilities 

liquidity 

R23 

Quick Assets/Net Sales 

liquidity 

R24 

Quick Assets/Total Assets 

liquidity 

R25 

Rate of Return to Common Stock 

profitability 

R26 

Retained Earnings/Total Assets 

profitability 

R27 

Return on Stock 

profitability 

R28 

Total Debt/Total Assets 

solidity 

R29 

Working Capital/Net sales 

liquidity 

R30 

Working Capital/Equity 

liquidity 

R31 

Working Capital/Total Assets 

liquidity 


prediction and their weights are obtained by applying multivariate discrimi¬ 
nate analysis. The final discriminate function obtained by Altman (say Za) is 
given by 

Za — 1-2 X\ 1.4 X 2 T 3.3 £3 0.6 X 4 -f- 0.999 x§. (2) 

The present condition of the firms may be assessed based on the Za values as 
follows: 

( Safe Zone, Za > 2.99; 

Za = < Grey Zone, 1.81 < Za < 2.99; 

[ Distress Zone, Za < 1-81. 

The companies in the safe zone may be considered financially healthy where 
as those in grey zone could go either way and if in the distress zone there is 
a greater risk that the company will go bankrupt within two years. Altman 
(2000) refined the weights of the Z-scores by revisiting the Zeta analysis and 
obtained the discriminate function as 


Z\j = 0.72 xi + 0.85 X 2 + 3.1 X 3 + 0.42 X 4 + x$. 


( 3 ) 
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In a linear model, the financial ratios Xi would influence Z -score in a lin¬ 
ear way. In the context of risk a change of a financial ratio by 1% may not 
have the same influence on the score (Tony et al. 2005; Atiya 2001; Baesens 
et al. 2003). Also, the studies in (Altman & Edith 2006; Aasen 2013) have 
found that financial ratios may be overstated due to accounting practices or 
manipulations. Therefore estimating nonlinear transformations for some of the 
independent variables would improve the bankruptcy predictions. The power 
transformations due to Box and Cox (Box & Cox 1964) are a popular method 
for nonlinear transformations to improve the symmetry and normality of the 
model fit. However, these transformations are proposed only on positive val¬ 
ues whereas the financial ratios could be negative. An alternative family of 
transformations is proposed in (Yeo & Johnson 2000) which may be applied 
even on negative values. 

In the present work we propose a nonlinear mapping Xi H > f(xi) of the 
form 



for transforming financial ratios (xt) before deriving the Z-score. The loglin- 
ear models make the differences between large values less important and those 
between small values more important. They are employed in (McLeay et al. 
2002; Ashton et al. 2004) to assume proportionate growth in accounting vari¬ 
ables that may be restricted to firm growth and are employed in (Senteney et 
al. 2006) for predicting the impending bankruptcy. 

We first transform the financial ratios utilizing the nonlinear function given 
in Equation 4 and then compute the new Z-score, Zm as 

Zm = Ai f(xi) + X 2 f (x 2 ) + A 3 /(x 3 ) + A 4 /( 2 : 4 ) + A 5 /(a^) +...+ A tf(xt), (5) 

in which Ai... A* denote the parameters of the financial ratios x±.. .Xt respec¬ 
tively. 

A more generalized form of Zm is given by 



( 6 ) 


in which Xk denotes the financial ratio for each k = 1 , 2 , ...,t and A k de¬ 
notes the weight of the Xk■ The weights A k are estimated using multivariate 
discriminate analysis (MDA). 

3 A New Indexing Measure for Credit Rating 


To predict bankruptcy of a firm it is required to identify the bounds on Z-score 
that can be estimated through empirical studies. These bounds are not compa¬ 
rable and they vary from business to business and also on countries’ economic 
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situations. Therefore, there is a need to standardize the Z-score to a distribu¬ 
tion before deriving useful indices for predicting the credit risk and bankruptcy 
of the firms. The present procedure generalizes the bankruptcy prediction that 
can be utilized by all the companies around the world. In (Hans et al. 2007) full 
credibility theory approach is used to estimate the parameters of the frequency 
(Poisson) and severity (Pareto) distributions for low frequency, high impact 
operational risk losses exceeding some threshold for each risk cell. An extreme 
value theory that provided fundamentals needed for the statistical modelling 
has been applied in stock market indices (Gilli & Kllezi 2006) to compute tail 
risk measures and extreme market events (Carvalhal & Mendes 2003). A new 
set of assessment models for long-term credit risk is studied in (Kubo & Sakai 
2011 ) which does not include stock prices and incorporates business cycles. 

The final objective of any risk model is to build the probability density 
function (PDF) of future losses in a portfolio. Renzo et al. (2006) developed 
simplest model using Bernoulli-distributed events and Poisson distribution. 
Probability distributions such as Poisson and Gamma have been employed to 
analyze aggregate loss distributions associated with operational risk (Degen 
2006; Dutta & Perry 2006; Embrechts et al. 2006). 

Creditrisk + models assume that the risk factors are independent gamma 
distributed random variables with mean 1 and variances a 1 (Matthias & 
Alexander 2012). Whereas, the P3 distribution with three parameters, loca¬ 
tion, shape and scale improves the goodness of fit to the data and can provide 
better estimates on the ratings. Moreover, Pearson family of distributions is 
employed in a wide range of applications such as financial time series mod¬ 
eling (Stavros 2014), distribution of stock returns (Pizzutilo 2002), flood risk 
modeling due to climate change (Miley et al. 2001). This distribution can fit 
a wide range of shapes with positive or negative skewness including a good 
approximation to the normal distribution. This motivated us to develop a 
new methodology for credit risk rating using three parameter P3 distribution. 
In addition, our methodology is universal and can apply to a wide range of 
distributions that are popularly being employed in credit risk applications. 

Presently, methods employed in credit risk applications use conventional 
moments such as mean, standard deviation, skewness and kurtosis to fit a 
distribution to the observations. These moments involve nonlinearities that 
are influenced by the presence of outliers and would result in over or under¬ 
estimation of the credit ratings. Therefore, we propose an approach based on 
the method of linearized moments popularly known as L-moments where in 
the parameters of the distribution can be expressed in a linear form. These 
L-moments can be computed using probability weighted moments (PWM) pre¬ 
sented in Section 3.1. The parameters of the P3 distribution (Hosking 1989) 
from L-moments can be computed using the procedure discussed in Section 3.2. 

In our methodology we first propose to fit a P3 distribution to the Zm score 
by computing L-moments and the parameters of the distribution. An index 
is then computed by measuring the deviations of the data using parameters 
of the P3 distribution. The ratings are then obtained by classifying the index 
into intervals ranging from highest safety (AAA) to high risk (CCC) based 
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on whether the value of the index is on positive or negative extreme of the 
distribution respectively. The details of the computations are presented in the 
Section 4. 


3.1 Probability Weighted Moments and L-moments 

The L-moments are analogous to conventional central moments, but can be 
estimated by linear combinations of order statistics. The L-moment estimates 
are found to be more robust compared to the conventional moments in the 
presence of outliers (Sankarasubramanian & Srinivasan 1999; Royston 1992; 
Ulrych et al. 2000). The L-moments are less sensitive to the effects of sam¬ 
pling variability, and are used to characterize a wide range of distributions 
than the conventional moments. Practically, they are less subject to bias in 
estimation and they approximate their asymptotic normal distribution more 
closely. The parameters estimated through L-moments are more accurate than 
the maximum likelihood and least square estimates. The L-moment estimates 
using PWM has been used in applications such as floods (Tai et al. 2012), 
drought (Eslamian et al. 2003) and financial risk (Maillet & Michel 2003). 

The probability weighted moments are defined in terms of the cumulative 
distribution function F(y ) (Greenwood et al. 1979) 

M P , r ,s = f F~ 1 (y) p F(y) r (l - F{y))'dF, (7) 

J o 

in which p, r, and s are positive integers, F^ 1 (y) denotes the inverse cumulative 
distribution function of the random variable Y . The term M P: r :S can now be 
used for describing the probability distribution. In a particular case where 
p = 1, and s = 0 the variable y becomes linear and the moment f3 r is defined 
as 

& = M 1|f ,o = / F~ 1 (y)F(y) r dF. (8) 

Jo 

The first three L-moments expressed as the linear combinations of the PWM 
as 


0 1 = A), (9) 

0 2 = 2/?i — /?0; 

03 = 6/3 2 - 6/3i + An 

in which 9\ known as L-mean, is a measure of central tendency and 02 known 
as L-standard deviation is a measure of dispersion. The ratios of L-moments 
are defined as 

r 2 =0 2 /0i, (10) 

T 3 = 63/62, 

in which 72 is termed as L-coefficient of variation and 73 is known as L-skewness 
and they are employed in estimating the parameters of the P3 distribution. 
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3.2 Pearson Type 3 (P3) Distribution 

In particular, the P3 probability density function g of the random variable E 
is defined as 

g(Z) = ^L[a(Z-c)y-'e-<*^ c \ ( 11 ) 

in which c, a and ry are location, scale and shape parameters of the distribution 
respectively. When the parameter a > 0, £ has positive skewness leading to 
c < £ < Too and when a < 0, £ has negative skewness leading to —oo < £ < 
c. Hence, c is a lower bound for positively skewed and an upper bound for 
negatively skewed P3 random variable E. 

The parameters c, a and ?y of P3 distribution are related to the L-moments 
as 

0.3333; 

r 3 < 1 


{ 1+0.2906 a 

5+0.1882 <5 2 + 0.04 42 <5 3 > 

0.36067 C—0.5967 ( 2 +0.2536 C 3 

1-2.78861t+2.56096C 2 -0.77045 C 3 ’ 


0 < r 3 < 
0.3333 < 


a = y/^e 2 e( r W- r ^+ 0 - 5 )), 
c = 6x - (a ?y), 

in which 5 = 3 7 rr| and £ = 1 — 73. 


( 12 ) 


4 Methodology 

The steps involved in the computation of the index are shown in Fig. 1. The 
dataset consisting of the financial ratios is first transformed into new variables 
using the nonlinear function proposed in Equation 6 . In the next step we con¬ 
vert the credit ratings given in the dataset into a binary variable henceforth 
known as bankruptcy index as follows. Define 


_ / 1 , V Rt e {B, BB, BBB, CCC }; . . 

\O,Vi? 0 e {A 7 AA,AAA}, 

in which R$ is the credit rating of the record (j) in the data set of m records, 
i.e., 4> takes values from l,...,m. As per CRISIL the credit ratings AAA 
denotes highest safety, AA denotes high safety, A denotes adequate safety, 
BBB denotes moderate safety, BB denotes moderate risk, B denotes high 
risk and CCC denotes very high risk. Clearly, from the Equation 13 one can 
infer that 6 ^ = 1 => bankruptcy or high risk category and = 0 =$■ non¬ 
bankruptcy or high safety category. 

Subsequently, the weights Ai,..., At of the transformed financial ratios are 
estimated using MDA with 6 ^ as dependent variable and f(x 1 ), /( x 2 ), ■ - ■, f(x t ) 
as independent variables. 
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bi — Ai /( Xi)i + X 2 f(x 2)1+, ■ ■ ■, +Xt f(xt) 1, 
b<j, = Ai f(x 1)0 + A2 f(x 2 )<f ,+, • • •, +A t /(x t )0, 


— ^1 /(®l)m + A 2 f(x 2 ) mH - ? • • • > + A( f{x t )m, (14) 

in which the variables f(x 1 ), f(x 2 ),..., f(xt ) denote the financial ratios after 
the application of the function / on the financial ratios * 1 , x 2 , ■ ■ ■, x t respec¬ 
tively as defined in Equation 4. The weights obtained from the Equation 14 
are substituted in Equation 5 to obtain the score Zm- 

The new Z-score Zm is then split into subsets Zm = ^ 2 ,{i 

..., Zfi where in j denotes the year of observation for the i th industry 
type. For each of these subsets the PWM are computed using the Equation 
8 . The L-moments di,d 2l d^, L-moment ratios t 2 and T 3 are computed using 
Equations 9, 10 respectively. The parameters c, 77 , a of P3 distribution are 
obtained from L-moments and L-moment ratio (73) by employing Equation 
12 . 



Fig. 1 Methodology for computing the new index measure 
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In the next step we standardize the dataset Zmi , with respect to the origin 
(parameter c) of the P3 distribution as Vij = ( Zm, ■ — c)/a. The new index 
Hij is then obtained as Hij = ({vij/ij) 033 + 1/(9 rf) — l) ( 977 )°'°. Based on 
the index the credit ratings are assigned to the dataset. The details of the 
procedure are presented in the Algorithm 1. 
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Algorithm 1 Algorithm for Computing Ratings based on the novel Credit 
Risk Index_ 

Require: 1. Data set S(m,n) where m and n denotes the number of rows and number of 
columns respectively. The attributes consist of financial ratios xi(m, 1), ..., xt(m , 1) 
and credit ratings R(m, 1) as per CRISIL. Clearly the number of columns in the 
dataset S is n = t + 1 . 

2. A set E(l,m) consisting of the type of industry type as in Ameya (2013) numbered 
from 1,..., 12 to which each row in the data set S belongs. 

3. A set Y(l,m) consisting of the year in which the ratings are observed for each row 
in the dataset S. 

Ensure: 1. Ratings W(l,m) for each row in the dataset S. 

Algorithm 

1. Transform the financial ratios given in the dataset S using the nonlinear transforma¬ 
tion function / given in Equation 4 to obtain a set with columns f(x i ),..., f(xt). 
Designate the set as D(m, t). The set D has m records and t attributes. 

2. Convert the credit ratings in column R in the dataset S to a bankruptcy index variable 
b using the transformation in Equation 13. We designate this set as 0. Obtain the 
set © = © U D. The set ©(m, n) has m rows and h = t + 1 columns. 

3. Obtain the weights Ai,..., A* i n Equation 14 by employing MDA with column b of 
the set © as dependent variable and the columns f(xi),...,f(xt) of the set © as 
independent variables. 

4. Obtain the new Z-score using the computations given in Equation 5 for all the records 
in the dataset © with the weights obtained in Step 3. Designate the set as 

Obtain the set Cl = EUYU@UZm- Clearly the set fi(m, ft) has m rows and h = n+3 
columns. 

5. For each type of industry i in column E of set fl 

(a) identify and collect all records in belonging to a particular type of industry 
E(i). The subset is designated as Q i.e Qm,n = m,n '■ wi £ E(i)}. Clearly 
Q C ft and rh < m. 

(b) obtain the PWM using the Zy column in set Q and then compute L-moments 
and their ratios using Equations 9, 10. 

(c) compute the parameters c, a, rj of P3 distribution using Equation 12. 

(d) for each year j in column Y of set Q 

i. compute the quantity Vij = (Zj — c)/a. 

ii. obtain the index Hij = {{vij/ 7 ^) 0 - 33 + 1 /( 977 ) — l) ( 977 ) 0 ' 5 . 

iii. obtain the row identification ( d ) corresponding to the i th industry and j th 
year i.e d = Index(i, j) would return the row number in dataset Q for which 
the rating is being computed. 

iv. compute the rating as defined by 


AAA, 

H itj > 2.0; 

AA, 

1.5 < Hij < 2.0; 

A, 

0 < Hij‘< 1.5; 

BBB, 

-1.0 <Hij < 0.0; 

BB, 

-1.5 < < -1.0; 

B, 

-2.0 < H-j < -1.5; 

CCC, 

H itj < -2.0, 


in which W denotes the set containing the new ratings, 

(e) initialize the set i.e Let Q = {0}- 

6. RETURN W. 

7. END. 
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4.1 Toy Example 

To illustrate our methodology, we utilize a toy dataset given in Table 2 with 
to = 10 records, financial ratios in columns x \,... ,xs (say t = 5), and the 
credit ratings in column R as per CRISIL. Clearly the total number of at¬ 
tributes n = t + l = 5 + l = 6. The dataset is designated as S(m, n). 


Table 2 Toy dataset 


Xi 

X 2 

X3 

X 4 

x 5 

R 

0.121 

0.263 

0.046 

1.219 

0.286 

BBB 

-0.046 

-0.164 

0.027 

0.218 

0.103 

B 

0.481 

0.696 

0.099 

3.969 

0.532 

AAA 

0.351 

0.238 

0.07 

1.023 

0.237 

BBB 

0.217 

0.326 

0.045 

2.522 

0.295 

AA 

0.105 

0.236 

0.053 

1.566 

0.216 

BBB 

0.078 

0.157 

0.041 

1.402 

0.335 

BBB 

0.189 

0.437 

0.059 

5.043 

0.452 

AAA 

0.043 

-0.047 

0.041 

0.287 

0.114 

B 

0.17 

0.702 

0.089 

23.002 

1.183 

AAA 


The dataset can as well be written as 

S = {(0.121,0.263,0.046,1.219,0.286, BBB), (-0.046, -0.164,0.027,0.218,0.103, B), (0.481,0.696,0.099,3.969,0.532, A 

(16) 

Clearly the cardinality of above set is |S| = 10. We now present the steps 
involved in Algorithm 1 for computing the credit ratings. We consider sets 
E(l, m) = {1,1,1,1,1,1,1,1,1,1} and Y(1,to) = {1,2,3,4,5,6,7,8,9,10} de¬ 
noting the type of industry and the year of observation respectively, for each 
of the records in S. Clearly the cardinality of the sets E and Y is equal to 
the number of rows in S i.e., |E| = |Y| = 10. A log-transformation function / 
given in Equation 4 is applied to each of the financial ratios (i.e for column X\ 
of record 1 we have f(x ij = /(0.121) = ln(0.121 + 1) = 0.114). We designate 
this set as D. Clearly the elements of this set are 

D = {(0.114,0.233,0.045,0.797,0.252), (-0.045, -0.152,0.027,0.197,0.098), (0.393,0.528,0.094,1.603,0.427), (0.301,0.2 

(17) 

We obtain the bankruptcy index b by applying Equation 13 to the column R 
in the Table 2 (i.e. record 1, Ri € BBB => bi = 1). We designate this set as 
0 and the elements of this set are {1,1, 0,1, 0,1,1, 0,1,0}. We then obtain a 
set 0 = DU0. The elements of this set are © are 

= {(0.114,0.233,0.045,0.797,0.252,1), (-0.045, -0.152,0.027,0.197,0.098,1), (0.393,0.528,0.094,1.603,0.427,0), (0.301 

(18) 

We obtain the set of equations between the bankruptcy index b and the 
financial ratios Xi, ... ,x§ as given in Equation 14. 
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b\ = 1 = 0.114 Ai + 0.233 A 2 + 0.045 A 3 + 0.797 A 4 + 0.252 A 5 , 
b 2 = 1 = -0.045 Ai - 0.152 A 2 + 0.027 A 3 + 0.197 A 4 + 0.098 A 5 , 

6 10 = 0 = 0.157 Ai + 0.532 A 2 + 0.085 A 3 + 0.178 A 4 + 0.781 A 5 . (19) 

By employing MDA on Equation 19 we obtain the weights as A 4 = 1.841, 

A 2 = -0.856, A 3 = -1.087, A 4 = 3.390, A 5 = -1.649. 

The score Zm is computed using Equation 5 as 

Z Ml = 1.841 x 0.114 - 0.233 x 0.856 - 0.045 x 1.087 
+0.797 x 3.390 - 0.252 x 1.649 = 2.249, 

Z Ml0 = 0.157 x 0.114 - 0.532 x 0.856 - 0.085 x 1.087 

+0.178 x 3.390 - 0.781 x 1.649 = 9.228. (20) 

The computed Zm scores obtained for the years j = 1,..., 10 using Equa¬ 
tion 20 is given as Z M = {2.249,0.525,4.900^2.335,3.914, 2.818, 2.464, 5.429,0.750,9.228}. 
We then obtain the dataset 0 = EUY U0U Zm as given in Table 3. 


Table 3 Data set obtained after applying function / in Equation 4 and bankruptcy index 
b in Equation 13 and new Z-score Zm using data in Table 2 


E 

Y 

f( x 1 ) 

/O 2 ) 

/O 3 ) 

f( x 4) 

/ 05) 

b 

Zm 

1 

1 

0.114 

0.233 

0.045 

0.797 

0.252 

1 

2.249 

1 

2 

-0.045 

-0.152 

0.027 

0.197 

0.098 

1 

0.525 

1 

3 

0.393 

0.528 

0.094 

1.603 

0.427 

0 

4.900 

1 

4 

0.301 

0.213 

0.068 

0.705 

0.213 

1 

2.335 

1 

5 

0.196 

0.282 

0.044 

1.259 

0.259 

0 

3.914 

1 

6 

0.100 

0.212 

0.052 

0.942 

0.196 

1 

2.818 

1 

7 

0.075 

0.146 

0.040 

0.876 

0.289 

1 

2.464 

1 

8 

0.173 

0.363 

0.057 

1.799 

0.373 

0 

5.429 

1 

9 

0.042 

-0.046 

0.040 

0.252 

0.108 

1 

0.750 

1 

10 

0.157 

0.532 

0.085 

3.178 

0.781 

0 

9.228 


The PWM are computed using Equation 8 to obtain the values /?o = 3.461, 
/3i = 2.449, and /3 2 = 1.939. 

The L-moments and their ratios are computed from /3o,/?i,/5 2 using Equa¬ 
tion 8, 10 as 

0 i = Po = 3.461, 

8 2 = 2 x 2.449 - 3.461 = 1.437, 

8 3 = 6 x 1.939 - 6 x 2.449 + 3.461 = 0.401, 
r 2 = 82/83 = 1.437/3.461 = 0.415, 
r 3 = 83/82 = 0.401/1.437 = 0.279. 


( 21 ) 
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The parameters of P3 distribution is obtained by substituting from L- 
moments obtained in Equation 21 in Equation 12. Since, T 3 < 0.333 we apply 


5 = 3 x 3.146 x (0.279) 2 = 0.7202, 

(1 + 0.2906 x 0.7202) 

V ~ (0.7202 + 0.1882 x (0.7202) 2 + 0.0442 x (0.7202) 3 ) 

(1 + 0.2093) 

= --- - - = 1 449 

(0.7202 + 0.0976 + 0.0165) ' ’ 

a = V3.1416 x 1.437 x e (T(i-449)-r(i.449+o.5)) 

= 1.7725 x 1.437 x e (-o.i2i4-(-o.o205)) 

= 2.5488 x e~° 1009 = 2.30 42, 

c = 3.461 - 1.449 x 2.3042 = 3.461 - 3.3398 = 0.121. (22) 

To compute the Index H for i = 1 and j = 1 we first compute 

t>i,i = (Zl hA -c)/a= (2.249 - 0.121)/2.3042 
= 2.1273/2.3042 = 0.9232, 

H hl = ((ui,i/1.449) 0 33 + 1/(9 x 1.449) - l) (9 x 1.449) 0 5 
= ((0.9232/1.449) 0 ' 33 + 1/(9 x 1.449) - l) x (9 x 1.449) 0 5 

= (0.8604 + 0.0767 - 1) x 3.6118 
= -0.0629 x 3.6118 = -0.2272. (23) 

The remaining indices of H for i = 1 and j = 2,..., 10 can be obtained 
using the steps in Equation 23 as Hi j 

= {-1.549,0.735, -0.186,0.433,0.028, -0.126,0.880, -1.265,1.711}. (24) 

To compute Wi we apply Equation 15 on H 11 = —0.2272 and find that it 
falls in the range —1.0 < —0.2272 < 0.0, hence the rating BBB is assigned to 
Wx. 

Similarly, the credit ratings for i = 1 and j = 2,..., 10 can be obtained 
using Equation 15 on values in Hij as W 2 i ... j io 

{B,A,BBB,A,A,BBB,A,BB,AA}. (25) 


5 Experiments and Results 

In this section we present the experiments conducted on the dataset and 
a comparison of our results with those of the earlier studies. A time series 
dataset Kubo & Sakai (2011) consisting of 3932 records with seven attributes 
as given in Table 4 is considered in our analysis. We have considered five fi¬ 
nancial ratios namely (i) working capital/total assets (WC_TA), (ii) retained 
earnings/total assets (RE_TA), (iii) earnings before interest and taxes/total 
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assets (EBIT_TA), (iv) market value of equity/book value of the total debit 
(MVE_BVTD) and (v) sales/total assets (S_TA) for the analysis. There are 
2392 bankrupt cases with credit ratings from BBB to CCC and 1540 non¬ 
bankrupt cases with credit ratings A to AAA. 


Table 4 Description of the dataset 


Sno 

Attribute 

Description 

Type 

1 

WC_TA 

working captial/total as¬ 
sets 

Real 

2 

RE TA 

retained earnings/total as¬ 
sets 

Real 

3 

EBIT_TA 

earnings before interest 
and taxes/total assets 

Real 

4 

MVE_BVTD 

market value of eq¬ 
uity/book value of total 
debit 

Real 

5 

S_TA 

sales/total assests 

Real 

6 

Industry 

1 to 12 

Categorical 

7 

Rating 

A, AA, AAA, B, BB, BBB, 
CCC 

Categorical 


The dataset consists of the credit ratings belonging to 12 different indus¬ 
tries with seven ratings ranging from highest safety (AAA) to very high risk 
(CCC). 


5.1 Results 

In this section we present the results of our study on the dataset. A comparison 
of skewness and the kurtosis of the original and transformed financial ratios 
are shown in Table 5. 


Table 5 Comparison of skewness and kurtosis of old and transformed financial ratios 


Financial 

ratio 

Skewness 

old 

Skewness 

transformed 

Kurtosis 

old 

Kurtosis 

transformed 

WC_TA 

-1.152 

-0.458 

17.944 

4.637 

RE.TA 

-2.476 

-1.591 

17.181 

6.462 

EBIT_TA 

-4.665 

-3.760 

74.310 

51.487 

MVE_BVTD 

12.992 

1.415 

269.574 

3.357 

S TA 

9.160 

2.129 

206.135 

12.598 


From Table 5 we can infer that the log transformation has reduced the 
skewness and kurtosis of the original variables, thereby improving the normal¬ 
ity of the financial ratios. 

We then estimate the weights of the loglinear model by performing MDA 
analysis between the bankruptcy index obtained by applying Equation 13 and 
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log transformed financial ratios. We have obtained the parameters given in 
the Equation 6 as Ai = 0.375, A 2 = 0.028, A 3 = —0.316, A 4 = 1.126, and 
A 5 = -0.236. 


Altman’s Z-score (Za) and the revised Z-score (Zjj) is computed as given 
in Equations 2 and 3. The Zm score of the transformed variables is computed 
and comparison of the descriptive statistics among the Zm , original Z-score 
and revised Z-score score is shown in Table 6 . 


Table 6 Descriptive Statistics among the new score (Zm), Altman’s Z-score (Za) and 
revised Z-score (Zu) 


Scoring 

Method 

Mean 

Standard 

deviation 

Quantile 

(25%) 

Median 

Quantile 

(75%) 

Zm 

1.007 

0.686 

0.529 

0.883 

1.374 

Za 

2.172 

3.060 

0.926 

1.626 

2.679 

Z U 

1.609 

2.178 

0.714 

1.207 

1.973 


We performed F-Test on the standard deviation of the Zm and other Z- 
score methods. The estimated F value ( F ca i ) is 30.664 and the tabulated F 
value (F tab ) is 39.863 at 0.01 significance. Since F cat is less than F tab we accept 
the null hypothesis that the two standard deviations are equal. The original 
Za and Zm are found to have good correlation with Spearman rho coefficient 
0.963 significant at 0.01 (2 tailed) level (see Table 7). 


Table 7 Correlation coefficients between Zm and Za 



Zm 

Z A 

Zm 

1.0 

0.963 

Za 

0.963 

1.0 


The estimates of the parameters of the P3 distribution obtained using 
the methodology described in Section 4 are shown in Table 8 for each of the 
industry types. 
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Table 8 P3 distribution parameter estimates using probability weighted moments 


Industry Type 

Location (e) 

Shape (/3) 

Scale (a) 

i 

-0.120815 

0.434659 

2.588683 

2 

-0.141720 

0.387980 

3.053867 

3 

-0.466324 

0.279429 

5.150730 

4 

-0.168323 

0.422125 

2.898292 

5 

-0.157205 

0.466022 

2.546405 

6 

-0.469409 

0.290401 

5.071536 

7 

-0.640034 

0.261942 

6.265737 

8 

-0.279976 

0.341867 

3.639809 

9 

-0.062825 

0.456268 

2.352211 

10 

-0.117848 

0.396729 

2.796561 

11 

-0.303738 

0.344976 

3.710743 

12 

-0.186726 

0.382682 

3.192083 


The credit ratings are computed after obtaining the standardized index 
of the data set using the P3 distribution parameters. The credit ratings are 
then converted to bankruptcy index using the Equation 13. A binary logistic 
regression, classification is carried out with the bankruptcy index as dependent 
variable and the Za, Zm or the Zjj scores as an independent variable as given 
in Equation 26 


A a = v a + v\ Za + Cl i 

A m = V m + v 2 Zm + ^2, (26) 

A u = V u + Z V + € 3 , 

in which A a , A m , A u denotes the bankruptcy index of the Z-scores Za , Zm 
and Zjj respectively. 

The parameters obtained from the three models are shown in Table 9. 


Table 9 Comparison of Z& , Zjj and the new score Zm for bankruptcy prediction using 
Wald statistics 


Parameter 

Za 

Zm 

Zu 

V a 

11.016(594.168) 

~ 

~ 

I'm 

- 

77.15(106.313) 

- 

IZu 

- 

- 

11.956(574.677) 

VI 

-5.423(572.703) 

- 

- 

"2 

- 

-78.534(106.569) 

- 

V3 

- 

- 

-7.993(555.379) 


The estimated coefficients of Za and Z M are both negative and statically 
significant at 0 . 01 % indicating that both measures are useful in predicting 
bankruptcy risk and lower the score, the higher the risk of bankruptcy. The 
coefficient of Zm is far lower than Za indicating the fact that the predictive 
power of Zm is far better than Altman’s Za score and revised Altman’s Zjj 


score. 
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A hold-out classification between Zm, Za , Zjj and bankruptcy index b is 
carried out using MDA and the prediction accuracy of the proposed method¬ 
ology is found to be 98.6% which is higher by 5% than any of the models 
proposed by Altman. This confirms that the proposed methodology is univer¬ 
sal and serves as a generalized tool that can improve the estimations of the 
existing methods/procedures in vogue and can predict the bankruptcy risk in 
an effective manner. 

Discriminate analysis on the data set generated using the new transfor¬ 
mation Zm has resulted in an accuracy of 93.7% in cross validated grouped 
cases correctly classified where as Altman’s Z-score Za has resulted only in an 
accuracy of 87.4%. 

MDA is carried out on Zm score and the credit ratings obtained from P3 
and Pareto distributions. The proposed method with P3 distribution resulted 
in an accuracy 92.2% whereas the model with the Pareto distribution has 
resulted in only 80% accuracy. 

To understand the sensitivity on the choice of thresholds for different rat¬ 
ings in predicting bankruptcy, we first construct a classification matrix or 
accuracy matrix (Table 10) based on the number of agreements and disagree¬ 
ments between the predicted group membership (estimated from the model) 
and the actual group membership (as present in the dataset) of bankruptcy 
for the thresholds given in Equation 15. 


Table 10 Classification matrix for the thresholds given in Equation 15 


Actual group membership 

Predicted group membership 

Bankrupt 

Non- B ankr up t 

Bankrupt 

Non-Bankrupt 

1966 (Ah) 
14 (M 2 ) 

426 (Mi) 

1526 (Ah) 


The actual group membership is equivalent to the a priori grouping and 
the predicted group refers to the cases wherein the proposed methodology 
attempts to classify them correctly. In the Table 10 N\, N 2 denotes the correct 
classifications (Hits) and Mi, M 2 denotes the misclassifications (Misses). N\ 
(1966) gives the number of cases of actual bankruptcy correctly classified as 
bankrupt by the proposed method. M\ (426) is the Type-I error that gives the 
number of cases wherein the actual group membership is bankrupt whereas 
the proposed model misclassified them as non-bankrupt. M 2 (14) is the Type- 
II error, that denotes the number of actual cases belonging to non-bankrupt 
group misclassified as bankrupt by the proposed model. N 2 (1526) are the 
number of cases wherein the proposed model correctly labels the actual cases 
as non-bankrupt. 

The accuracy of the proposed methodology is computed as (Ni+N 2 )/(Ni + 
N 2 + M 1 +M 2 ) = (1966+1526)/(1966+1526+14+426) = 3492/3932 = 0.88 = 
88 %. The Type-I error is the ratio of misclassified cases of actual bankrupt 
cases declared as non-bankrupt by the model with total bankrupt cases i.e 
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Type-I = M l /{N 1 + M{) = 426/(1966 + 426) = 426/2392 = 0.177 = 17.7%. 
The Type-II error is the ratio of misclassified cases of actual non-bankrupty 
cases declared as bankrupt by the model with total non-bankrupty cases i.e 
Type-II=M 2 /(A 2 + M 2 )=14/(14 + 1526) = 14/1540 = 0.009 = 0.9%. 

The proposed method with thresholds as given in Equation 15 is accurate 
in classifying 88.8% of total samples with Type I error to be only 17% while the 
Type II error was even better at 0.9%. Therefore, there is a positive upward 
bias which can be addressed by adjusting the thresholds between the credit 
ratings A and BBB as the boundaries fall in the grey zone. Keeping the other 
thresholds unchanged, we updated the thresholds of BBB as —1.0 < H, j < 
0.25 and A as 0.25 < H,j < 1.5 from the classification table, we obtained the 
Type I error as 4% and Type II error as 5% with overall accuracy of 95%. To 
see if the sensitivity be further improved we updated the thresholds of BBB 
as —1.0 < Hi j < 0.5 and A as 0.5 < Hjj < 1.5 keeping the others unchanged. 
We found from the classification table the Type I error as 0.16% whereas Type 
II error has increased to 21.7% with overall accuracy of 91.4%. Therefore, the 
choice of thresholds for transition from bankruptcy to non-bankruptcy should 
be chosen with caution so that both Type I and Type II errors are at minimum. 

Even though the samples are disproportionate the Algorithm 1 has out¬ 
performed the accuracies obtained using the Altman’s Z-score methods. 


6 Conclusions and Discussion 

A new nonlinear transformation procedure for building a log linear model for 
computing the Z-score is proposed. Based on the new Z-score (Z M ) a new 
indexing measure is proposed by fitting the data to a P3 distribution and 
then obtaining the deviations of the given dataset from the standard nor¬ 
mal using an equi-probability transformation. The multivariate discriminate 
analysis (MDA) for predicting the bankruptcy index has shown that the pro¬ 
posed methodology has given highest accuracy of 98.5% which is higher by 
5% as compared with Altman’s Z-score. The classification accuracies of the 
transformed financial ratios in predicting the bankruptcy is around 93.7% as 
compared to 87.4% obtained by using the factors of Altman’s procedure. The 
accuracies of the proposed method with P3 distribution was 92.2% where as 
a model with Pareto distribution resulted in an accuracy of 80%. Though the 
methodology is universal and serves as a generalized tool, there is an immense 
need to validate with global datasets. Also, mutual interference among finan¬ 
cial ratios is an important aspect that requires further investigation. We defer 
our ongoing work in this direction to a subsequent exposition. 
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